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1  Project  Objectives  and  Status 

1.1  Review  of  Program  Objectives 

The  base  effort  was  intended  to  address  the  deficiency  in  the  understanding  of  the  effect  that 
near-nozzle  and  inner-nozzle  flow  conditions  have  on  jet  noise  radiation.  There  were  three 
primary  objectives: 

1.  high-fidelity  characterization  of  heated,  over-expanded  supersonic  jets, 

2.  source  identification  through  development  of  advanced  analytical  diagnostics,  and 

3.  enhanced  computational  modeling  of  hot  supersonic  jets. 

1.2  Project  Status 

This  reporting  period  marks  the  end  of  the  base  effort.  Figure  1  shows  the  project  status. 
The  following  list  summarizes  completion/incompletion  of  specific  WBS  entries. 

WBS  1.4  The  methodology  for  computation  of  the  relevant  terms  was  completed  using 
data  generated  from  the  LES  (large-eddy  simulation)  in  WBS  2.3.  These  results  were 
presented  at  the  2013  AIAA/CEAS  Aeroacoustics  conference  as  papers  2013-2081  and 
2013-2193. 
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WBS  2.1.4  The  temperature  profile  measurements  were  not  completed.  The  first  attempt 
ended  due  to  probe  failure.  Focus  on  the  MHz  rate  measurements  in  WBS  2.2.2  took 
presedence  and  we  did  not  revisit  these  measurements. 

WBS  2.2.2  The  campaign  was  completed,  but  the  measurements  were  less  successful  than 
anticipated.  Details  are  discussed  in  Section  2.2. 

WBS  3.0  The  main  deficiencies  with  the  experimental  measurements  were  related  to  the 
Cordin  camera:  misalignment  of  the  individual  MCP  (multi-channel  plate)  units  and 
resonant  vibration  of  the  pellicle  beam  splitters.  Neither  of  these  problems  could  be 
rectified  in  the  lab.  As  a  result,  we  could  not  complete  the  related  data  analysis. 
Efforts  to  extract  at  least  some  velocity  data  are  being  conducted  by  Dr.  Thurow’s 
graduate  student,  Bryan  Brock,  and  will  be  presented  at  the  upcoming  AIAA  SciTech 
Conference. 

WBS  5.0  Plenoptic  PIV  data  was  acquired  to  attempt  a  different  method  of  flow  field 
measurement.  Results  and  suggestions  for  future  work  are  discussed  in  Section  2.2. 
WBS  6.0  The  phase-array  methodology  was  completely  implemented  and  results  were  pre¬ 
sented  at  the  peer  review. 

WBS  8.0  This  report  completes  the  quarterly  reporting  in  the  base  period. 

WBS  9.0  A  final  report  is  not  required  at  this  time  as  a  no-cost  extension  has  been  granted 
in  connection  with  the  project  extension. 

WBS  10.2  The  year  2  peer  review  was  completed.  Reviewer  comments  led  to  a  proposed 
reduction  in  scope  for  the  year  3  project  extension. 

1.3  Extension  Proposal 

Given  the  lack  of  progress  on  the  time-resolved  flow-field  measurements,  we  were  encouraged 
to  propose  a  modified  work  plan  for  a  contract  extension  that  focused  on  the  near-field 
propagation  of  jet  noise.  The  extension  work  will  address  two  question. 

1.  Do  the  shock-like  structures  form  inside  the  jet  flow  or  are  they  the  result  of  rapid 
steepening  immediately  outside  of  the  jet? 

2.  What  is  the  correct  decay  law  for  the  acoustic  signature? 

The  extension  proposal  was  submitted  on  16  September  2013  with  the  project  plan  shown 
in  Figure  2 

2  Activity  for  Current  Reporting  Period 

2.1  Jet  Acoustics  and  Crackle 

The  progress  on  the  non-linear  propagation  work  was  compiled  and  summarized  in  Woutijn 
Baars’s  dissertation  under  the  direction  of  Dr.  Tinney.  The  summary  is  provided  in  Ap¬ 
pendix  A.  The  summary  also  provide  suggestions  for  future  work  that  are  relevant  to  our 
efforts  moving  forward  in  the  extension  period  of  this  program. 


Page  3 

Quarterly  Progress  Report  No.  8  -  N00014-1 1-1-0752 

Approved  for  Public  Release  Distribution  is  Unlimited 


Active  Control  of  Noise  from  Hot  Supersonic  Jets  |  Quarterly  Rpt.  8 


Figure  1.  Project  chart  showing  WBS  items  and  current  completion  status  as  of  the  end  of 
the  reporting  period. 
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Extension  Work  Period 


1  Oct  2013  -30  September  2014 


WBS  1.0  Determine  Amplitude  and  Decay  Rate  of 
Acoustic  Waveform 


WBS  1. 1  Separating  hydrody  namic  and  acoustic  com- 
ponents  of  the  wave  form . 


WBS  1.2  Conduct  measurements  on  a  Mach  3  shock- 
free  jet  using  1  /8-inch  microphones. 


WBS  1.3  Develop  Go l  * dburg  numbers  for  jet  noise 
problems.   


WBS  1.4  Demonstrate  and  validate  the  decay >  law 
model  with  acoustic  data  from  hot,  shock-containing 
jets. 


WBS  2.0  Prepare  Quarterly  Reports 
WBS  3.0  Present  Results  at  Annual  Review  Meeting 
WBS  4.0  Prepare  Final  Report 


Q2 

Q3 

Figure  2.  Extension  work  plan. 


2.2  High-Fidelity  Flow  Field  Measurements 

2.2.1  Time- Resolved  PIV:  Lessons  Learned 

As  indicated  in  the  review  of  the  base  effort  and  evidenced  in  our  annual  reports,  progress 
with  time-resolved  PIV  lms  been  slow  and  frustrating.  In  the  first  year,  a  combination  of 
poor  particle  seeding  and  poor  image  quality  stifled  progress.  In  the  second  year,  particle 
seeding  issues  were  overcome;  however,  the  image  quality  and  slow  acquisition  rate  of  the 
Cordin  camera  severely  hampered  progress  and  limited  acquisition  to  a  single  set  of  data 
(over  300  image  sequences  were  acquired  and  are  currently  being  processed).  As  such,  the 
vast  majority  of  the  data  processing  effort  has  been  the  development  of  image  processing 
algorithms  to  overcome  these  issues  as  opposed  to  exploiting  the  unique  information  afforded 
by  this  system.  Unfortunately,  these  issues  have  prevented  us  from  proceeding  on  to  more 
interesting  aspects  such  as  asynchronous  triggering. 

2.2.2  Plenoptic  PIV:  A  Different  Experience  than  Time- Resolved  PIV 

In  the  2nd  year  of  the  base  effort  we  conducted  a  simple  proof-of-concept  experiment  in 
the  NCPA  facility  with  a  plenoptic  camera  designed  and  built  by  co-PI  Thurow  to  test 
the  feasibility  of  performing  3D  velocity  measurements  in  a  heated  supersonic  jet.  In  stark 
contrast  to  the  experience  with  time-resolved  PIV,  the  experience  with  the  plenoptic  camera 
went  very  smoothly  as  we  were  able  to  acquire  over  300  image  pairs  in  1  day  of  testing 
without  any  major  issues.  This  included  both  the  time  to  set-up  the  experiments  as  well  as 
acquire  the  data.  As  will  be  detailed  below,  this  simple  proof-of-concept  experiment  yielded 
good  quality  data  that  reveals  details  about  the  large-scale  3D  structure  of  the  jet  flow. 
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Figure  3.  Illustration  of  ability  to  computational  refocus  an  image  after  it  has  been  acquired. 

2.2.3  Plenoptic  PIV :  Background  and  Sample  Results 

As  plenoptic  PIV  was  not  part  of  the  originally  proposed  effort,  a  brief  background  of  the 
technique  is  provided  here.  More  details  are  available  in  [1,  2,  3,  4,  5].  A  plenoptic  camera 
is  distinguished  from  a  conventional  camera  in  its  unique  ability  to  record  the  4D  light- 
field  of  a  scene.  The  term  light-field  describes  the  complete  distribution  of  light  rays  in 
space  leading  to  a  4D  function,  sometimes  termed  the  plenoptic  function,  wherein  each 
light  ray  is  characterized  by  its  position  (x,y)  on  a  plane  and  its  propagation  angle  (#,0). 
Standard  imaging  captures  only  two  dimensions  of  the  4D  light-field  as  angular  information  is 
integrated  at  the  sensor’s  surface,  leading  to  familiar  imaging  phenomena  such  as  de-focus, 
depth-of-field  and  image  blur.  A  plenoptic  camera,  on  the  other  hand,  uses  a  microlens 
array  placed  near  a  conventional  image  sensor  to  measure  both  the  position  and  angle  of 
light  rays  captured  by  the  camera  lens,  thus  forming  a  4D  measurement  of  the  light-field. 
Co-PI  Thurow  has  recently  designed  and  assembled  a  plenoptic  camera  in  his  laboratory 
with  its  key  attributes  being  its  compact  size,  simple  design,  ease  of  operation  and  ability 
to  collect  light-field  data  in  a  single  snapshot.  Figure  3  illustrates  one  of  the  many  uses  of 
light-field  data:  the  ability  to  computationally  refocus  an  image  after  acquisition.  These 
images  were  computed  from  a  single  snapshot  using  an  in-house  ray-tracing  algorithm  to 
virtually  propagate  the  light  rays  forward  and  backward  in  space,  and  integrate  the  angular 
information  in  a  manner  similar  to  a  conventional  camera.  This  refocusing  feature  is  the 
most  noteworthy  aspect  of  the  consumer  camera  developed  by  Lytro  Inc.  with  variations  of 
this  technology  also  available  from  Raytrix.  Another  common  use  of  the  data  is  to  generate 
views  obtained  with  different  perspectives.  With  our  camera,  we  are  able  to  generate  over 
100  unique  views  from  different  angles  from  a  single  snapshot  of  data. 

The  emphasis  of  Co-PI  Thurow’s  work  over  the  last  three  years  has  been  on  adaptation 
of  this  technology  for  3D  velocity  measurements.  This  has  included  development  of  both  the 
hardware  and  software  necessary  to  realize  3D  velocity  measurements  in  a  turbulent  flow. 
The  first  experimental  demonstration  of  the  camera  focused  on  the  collection  of  plenoptic 
PIV  image  data  in  both  ZPG  (zero  pressure  gradient)  and  APG  adverse  pressure  gradi¬ 
ent  turbulent  boundary  layers  (Ree  =  4, 500  and  7,  200,  respectively).  To  determine  the  full 
3D/3C  velocity  field,  a  pair  of  3D  pmrticle  fields  were  reconstructed  from  individual  snapshots 
of  the  light-field  data  using  limited-angle  tomography  algorithms  (the  MART  algorithm  used 
in  tomo-PIV).  Conventional  3D  PIV  cross-correlation  algorithms  were  then  applied  to  the 
particle  volumes  to  determine  the  particle  displacement  and  velocity.  Examples  of  other 
flows  in  which  the  plenoptic  camera’s  capabilities  have  been  demonstrated  are  the  aforemen¬ 
tioned  turbulent  boundary  layers,  visualization  of  the  counter-rotating  vortices  generated  by 
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Figure  4.  Sample  3D  velocity  field  obtained  in  a  2"  diameter  supersonic  jet  seeded  with 
alumina  particles.  Y-axis  is  stretched  by  a  factor  of  5  to  illustrate  different  cross-sections  of 
the  jet.  Color  corresponds  to  streamwise,  uy ,  component  of  velocity. 


a  Wheeler  doublet  vortex  tab  in  a  turbulent  boundary  layer  and  the  supersonic  jet  results 
discussed  here. 

Figure  4  shows  an  example  of  the  data  recently  acquired  in  the  NCPA  demonstration. 
It  is  an  instantaneous  3D/3C  ( 3-dimension al/3-component)  velocity  field  obtained  in  the 
heated,  supersonic  jet  at  the  center  of  this  study.  The  jet  Mach  number  was  1.74,  the  total 
temperature  ratio  was  3.37  and  the  jet  diameter  50.8  mm.  The  dimensions  of  the  volume 
are  approximately  100  x  91  x  66  mm  centered  approximately  1.5  jet  diameters  downstream 
from  the  nozzle  exit.  The  volume  was  illuminated  with  a  200  mJ/pulse  PIV  laser  and  the 
camera  was  mounted  approximately  18”  away  from  the  jet.  Particle  seeding  was  obtained  in 
the  same  fashion  as  previous  PIV  experiments  using  alumina  particles  suspended  in  ethanol 
and  injected  into  the  jet’s  settling  chamber.  All  calculations  were  performed  on  a  standard 
desktop  computer. 

A  selection  of  image  pairs  from  the  November  2012  effort  were  analyzed.  A  subset  of 
24  image  pairs  were  processed  and  used  to  determine  mean  and  rms  quantities.  Figure  5(a) 
shows  the  average  velocity  field  calculated  from  this  limited  set  of  data.  A  key  feature 
of  the  average  velocity  field  is  the  presence  of  a  low*  velocity  region  in  the  center  of  the 
jet  resulting  from  the  slipline  where  shocks  converge.  This  matches  probe  measurements 
and  computational  simulations.  Other  features,  such  as  the  relatively  thin  shear  layer  and 
magnitude  of  the  velocity  is  also  consistent  with  expectations.  Figure  5(b)  shows  the  RMS 
fluctuations  of  the  velocity  field.  The  jet  core  and  ambient  both  showr  lowr  levels  of  fluctuations 
whereas  the  shear  layer  shows  the  highest  levels  of  fluctuations.  The  key  point  is  that 
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(a)  Mean  Velocity  Field  (b)  RMS  Velocity  Fluctuations 

Figure  5.  Average  and  RMS  velocities  calculated  using  24  image  pairs.  Five  cross-sections 
are  shown  with  the  lower  velocity  jet  core  apparent  near  the  center  of  the  jet.  Velocity  scale 
ranges  from  200  (blue)  to  800  (red)  m/s. 


plenoptic  PIV  accurately  shows  the  three  dimensional  location  where  velocity  fluctuations 
are  known  to  occur. 

2.2.4  Plenoptic  PIV:  Comment  on  Occlusions 

A  concern  was  expressed  about  how  the  plenoptic  PIV  technique  can  overcome  occlusions, 
presumably  due  to  the  notion  that  particles  wall  obstruct  the  view  of  other  particles.  We 
do  not  believe  this  to  be  a  problem  for  the  following  reasons.  In  principle,  for  very  low 
image  densities  where  particles  can  be  explicitly  identified,  the  3D  location  of  a  particle 
can  be  determined  using  only  two  images  acquired  from  different  viewing  angles,  such  as 
achieved  commonly  in  particle  tracking  velocimetry  experiments.  For  higher  image  densities 
where  individual  particles  cannot  be  explicitly  identified,  tomo-PIV  has  recently  emerged 
as  a  viable  technique  for  3D  flow  measurements.  I11  tomo-PIV,  tomographic  algorithms  are 
used  to  reconstruct  a  volume  of  particles  based  on  2D  projections  (e.g.  images)  of  the  volume 
onto  cameras  located  at  different  angular  positions  around  the  volume.  I11  this  configuration 
occlusions  are  still  present  for  some  of  the  viewTs;  however,  the  additional  views  provide 
enough  additional  information  for  the  tomographic  algorithms  to  accurately  reconstruct  the 
location  of  particles  within  the  volume.  The  effect  of  occlusions  is  a  reconstruction  artifact 
known  as  ‘ghost  particles’  that  degrades  the  quality  of  the  reconstruction.  It  was  found 
that  the  influence  of  ghost  particles  decreases  significantly  as  one  increases  the  number  of 
cameras  in  the  system,  with  four  cameras  considered  reasonable  for  many  applications.  We 
do  not  believe  plenoptic  PIV  is  affected  by  this  problem  because,  in  the  context  of  the 
above  description,  the  plenoptic  camera  produces  the  virtual  equivalent  of  over  100  distinct 
projections  of  the  volume.  We  have  tested  this  using  a  synthetic  image  generator  developed 
specifically  for  particle  imaging  with  a  plenoptic  camera.  Figure  6  shows  the  reconstruction 
quality,  Q<  as  a  function  of  the  particle  number  density.  In  this  case,  the  number  of  pixels 
is  taken  as  the  total  number  of  active  pixels  on  the  camera  which  is  equal  to  11  x  106.  As 
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Figure  6.  Reconstruction  quality  as  a  function  of  number  of  particles.  Particles  per  pixel  are 
calculated  based  on  the  total  number  of  active  pixels  on  the  camera  (11  x  106  pixels). 


can  be  seen,  the  quality  of  the  reconstruction  is  relatively  flat  over  a  large  range  of  particle 
number  densities  and  exhibits  a  graduate  decrease  for  large  particle  number  densities.  As 
such,  we  believe  plenoptic  PIV  is  fairly  robust  and  resistant,  particularly  in  comparison  to 
techniques  such  as  tomo-PIV. 

Along  the  same  lines,  plenoptic  PIV  affords  additional  advantages  over  tomo-PIV  that 
make  it  well  suited  for  this  study.  Typical  tomo-PIV  systems  are  limited  to  relatively  thin 
(order  10  mm  thick)  volumes  partially  to  limit  the  particle  number  density  of  the  projections, 
but  also  to  satisfy  the  requirement  that  the  full  volume  remain  within  the  depth-of-field  of 
each  camera.  To  achieve  this,  each  camera  lens  is  typically  adjusted  to  a  small  aperture 
setting  to  increase  the  depth-of-field.  This  significantly  limits  the  amount  of  light  collected 
by  the  camera,  further  restricting  the  size  of  the  volume  due  to  the  need  for  very  high  energy 
laser  pulse  illumination.  Plenoptic  PIV,  on  the  other  hand,  avoids  these  complications  as 
it  does  not  require  the  volume  to  be  “in-focus”  in  a  conventional  sense.  Rather,  plenoptic 
PIV  actually  works  better  with  a  wide  open  aperture  due  to  the  increased  angular  content 
sampled  by  the  camera.  This  was  illustrated  by  the  November  2012  NCPA  acquisition  under 
the  base  effort  where  a  conventional  PIV  laser  provided  sufficient  illumination.  Combined 
with  the  mitigation  of  ghost  particles,  plenoptic  PIV  is  well  suited  for  3D  measurements  of 
‘thick’  volumes. 

[1]  Fahringer,  T.  &  Thurow,  B.  (2012)  “Tomographic  Reconstruction  of  a  3-D  Flow  Field 
Using  a  Plenoptic  Camera,”  in  42nd  AIAA  Fluid  Dynamics  Conference  (New  Orleans,  LA). 

[2]  Fahringer,  T.  &  Thurow,  B.  (2013)  “The  Effect  of  Grid  Resolution  on  the  Accuracy  of 
Tomographic  Reconstruction  Using  a  Plenoptic  Camera,”  in  51st  AIAA  Aerospace  Sciences 
Meeting. 
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Imaging,  Master’s  thesis,  Auburn  University. 
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3  Technical/Cost  Status  &  Problem  Areas 

The  application  of  MHz  PIV  to  the  high-noise  environment  near  the  jet  proved  a  technical 

challenge  that  could  not  overcome  in  the  present  work.  As  a  result,  the  extension  work  plan 

was  reduced  to  focus  on  the  propagation  of  the  high  amplitude  jet  noise  acoustics. 

4  Publications,  Meetings,  and/or  Travel 

4.1  Running  List  of  Publications  Produced 

[1]  Baars,  W.J.,  Tinney,  C.E.,  Murray,  N.E..  Jansen,  B.J.  &  Panickar,  P.  (2011)  “The 
Effect  of  Heat  on  Turbulent  Mixing  Noise  in  Supersonic  Jets,”  m  4 9th  AIAA  Aerospace  Sciences 
Meeting ,  Paper  2011-1029  (Orlando,  FL). 

[2]  Baars,  W.J.,  Tinney,  C.E.  &  Wochner,  M.S.  (2012)  “Nonlinear  Noise  Propagation  from 
a  Fully  Expanded  Mach  3  Jet,”  in  50th  AIAA  Aerospace  Science  Meeting ,  Paper  2012-1177 
(Nashville,  TN). 

[3]  Murray,  N.E.,  Lyons,  G.W.,  Tinney,  C.E..  Donald,  B..  Baars,  W.J.,  Thurow,  B.S., 
Haynes,  R.H.  &  Panickar,  P.  (2012)  “A  Laboratory  Framework  for  Synchronous  Near/Far- 
Field  Acoustics  and  MHz  PIV  in  High-Temperature,  Shock- Containing,  Jets,”  in  Proceedings 
of  the  Internoise  2012/ASME  NCAD  Meeting ,  ASME/NCAD-1270  [invited]  (New  York  City. 
NY). 

[4]  Baars,  W.J.  &  Tinney,  C.E.  (2012)  “Scaling  Model  for  Nonlinear  Supersonic  Jet  Noise,” 
in  Bulletin  of  the  American  Physical  Society ,  57:17.  Abstract  D24:8  (San  Diego,  CA). 

[5]  Panickar,  P.,  Erwin,  J.,  SInha,  N.,  Murray,  N.E.  &  Lyons,  G.W.  (2013)  “Localization 
of  Acoustic  Sources  in  Shock- Containing  Jet  Flows  Using  Phased  Array  Measurements,”  in 
51st  AIAA  Aerospace  Sciences  Meeting ,  2013-0613  (Grapevine.  TX). 

[6]  Haynes,  R.H.,  Brock,  B.A.  &  Thurow,  B.S.  (2013)  “Application  of  MHz  Frame  Rate, 
High  Dynamic  Range  PIV  to  a  High- Temperature,  Shock-Containing  Jet,”  in  51st  AIAA 
Aerospace  Sciences  Meeting ,  2013-0774  (Grapevine,  TX). 

[7]  Fievet,  R.,  Tinney,  C.E.,  Murray,  N.E.,  Lyons,  G.W.  &  Panickar,  P.  (2013)  “Acous¬ 
tic  Source  Indicators  using  LES  in  a  Fully  Expanded  and  Heated  Supersonic  Jet,”  in  19th 
AIAA/CEAS  Aeroacoustics  Conference ,  2013-2193  (Berlin,  Germany). 
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[8]  Baars,  W.J.  &  TlNNEY,  C.E.  (2013)  “Quantifying  Crackle-Inducing  Acoustic  Shock- 
Structures  Emitted  by  a  Fully-Expanded  Mach  3  Jet,”  in  19th  AIAA/CEAS  Aeroacoustics 
Conference,  Paper  2013-2081  (Berlin,  Germany). 

[9]  Baars,  W.J.  &  TlNNEY,  C.E.  (2013)  **A  Temporal  and  Spectral  Quantification  of  the  Crackle 
Component  in  Supersonic  Jet  Noise,”  in  2nd  Symposium  on  Fluid-Structure-Sound  Interactions 
and  Control ,  p.  170  (Shatin,  Hong  Kong). 

[10]  Baars,  W.J.,  Tinney,  C.E.  &  Wochner,  M.S.  “Nonlinear  Distortion  of  Acoustic  Wave¬ 
forms  from  High-Speed  Jets,”  J.  of  Fluid  Mechanics  [in  review] . 

[11]  Baars,  W.J.  (2013)  Acoustics  from  High-speed  Jets  with  Crackle ,  Ph.D.  Dissertation,  The 
University  of  Texas  at  Austin. 

5  Planned  Activities  fur  Next  Reporting  Period 

•  Initiate  work  under  the  project  extension  and  complete  all  related  subcontract  modi¬ 
fications. 

•  Extract  a  data  set  from  the  LES  with  more  samples  and  a  larger  spatial  domain  than 
used  previously  for  analysis  under  the  project  extension. 

•  Conduct  measurements  on  a  Mach  3  shock-free  jet  using  1/8-inch  microphones. 
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''Summary  and  Future  Work,”  from  Acoustics  from  High-Speed  Jets 
with  Crackle  by  Woutijn  Baars 
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Chapter  7 

Summary  and  future  work 

The  summary  is  condensed  into  two  parts.  At  first,  the  research  on  cu¬ 
mulative  nonlinear  waveform  distortions  in  the  acoustic  field  of  high-speed  jets 
is  reviewed  (§  7.1).  Secondly,  the  investigation  of  crackle-inducing  signatures 
is  summarized  in  §  7.2.  Suggestions  for  future  work  are  itemized  throughout 
the  summary. 

7.1  Cumulative  nonlinear  waveform  distortions 

The  acoustic  field  of  an  unheated  and  perfectly  expanded  Mach  3  ex¬ 
perimental  jet  was  examined  to  understand  the  degree  of  both  local  and  cu¬ 
mulative  nonlinear  distortions  of  waveforms  emitted  by  high-speed  jets.  The 
source  mechanisms  of  interest  are  those  that  are  known  to  produce  Mach  waves, 
which  become  pronounced  as  the  convective  speed  of  the  large-scale  events  in 
the  shear  layer  move  at  supersonic  speeds  relative  to  the  ambient  surround¬ 
ings.  Most  intense  noise  is  emitted  from  a  region  bounded  by  the  potential 
core  and  the  point  of  sonic  velocity  ou  the  jet’s  centerline.  The  resultant  sound 
field  forms  a  highly  directive  pattern  with  a  ridge  that  follows  along  the  Mach 
cone  half  angle  of  the  jet. 
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A  model  for  predicting  the  presence  of  cumulative  nonlinear  effects, 
distorting  acoustic  waveforms  produced  by  jet  flows,  was  then  proposed.  The 
model  considers  the  broadband  sound  along  the  peak  noise  path  and  is  based 
on  an  assessment  of  the  effective  GoPdberg  number:  a  ratio  of  shock  formation 
distance  to  acoustic  absorption  length.  Relevant  scaling  parameters  are  shown 
to  comprise  Mach  number  (gas  dynamic  or  acoustic),  temperature  ratio,  and 
Strouhal  number  corresponding  to  the  peak  frequency  of  the  broadband  noise. 
The  model  assumes  that  the  noise  is  emitted  from  a  location  that  scales  with 
nozzle  diameter.  Model  predictions  are  justified  by  several  laboratory-  and  full- 
scale  data  sets.  Two  different  scaling  scenarios  are  presented  for  the  practising 
scientist  to  choose  from.  The  first  of  these  allows  shock  formation  distance  to 
be  preserved  between  laboratory  ami  full-scale  conditions  and  is  based  on  a 
scaling  of  the  shock  formation  distance  to  nozzle  diameter.  The  second  scenario 
allows  one  to  preserve  the  wave  steepening  process  between  laboratory  and  full- 
scale  conditions  based  on  a  match  of  the  effective  Gol  dberg  number.  Both 
methods  provide  promising  results  upon  application  of  the  laboratory  and  full- 
scale  experimental  data  sets.  There  are,  of  course,  opportunities  to  improve 
the  model;  these  are  proposed  in  itemized  form. 

■  A  more  thorough  assessment  of  how  jet  exit  conditions  affect  scaling 
parameters  has  to  be  pursued.  Where  Strouhal  number  is  concerned, 
it  is  still  unclear  how  the  peak  frequency  along  the  peak  noise  propa¬ 
gation  path  depends  on  jet  exit  conditions  such  as  temperature  ratio 
or  Mach  number.  It  is  also  assumed  here  that  the  location  from  where 
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the  waveform  is  emitted  scales  with  nozzle  diameter  only.  This  assump¬ 
tion  disregards  the  dependence  of  the  emission-point  location  on  Mach 
number,  Reynolds  number,  temperature  ratio,  or  even  the  operating 
state  of  the  nozzle  (overexpanded  or  underexpanded).  Deficiencies  in 
the  understanding  of  these  dependencies  would  be  alleviated  by  accu¬ 
rate  parametrization  of  the  effects  of  Strouhal  number,  Mach  number, 
temperature  ratio,  and  Reynolds  number  on  the  emission-point  location 
of  the  broadband  noise  that  is  an  input  to  the  model.  Taking  the  ef¬ 
fect  into  account  of  jets  operating  at  overexpanded  or  underexpanded 
conditions  are  of  great  practical  significance  to  this.  For  off-design  con¬ 
ditions,  broadband  shock  noise  is  saturated  by  turbulent  mixing  noise 
along  shallow  angles  to  the  jet  axis  (within  the  Mach  cone  of  the  jet)  so 
distortions  to  the  acoustic  waveform  are  still  anticipated  so  long  as  the 
effective  Gohdberg  number  is  shown  to  be  sufficiently  large. 

■  It  is  unclear  to  what  effect  a  pre-steepened  acoustic  waveform  has  on 
shock  formation  distance.  Waveform  steepening  occurs  when  the  wave 
amplitude  is  large  enough  to  overcome  viscous  absorption.  This  requires 
a  large  source  intensity  and  hence,  a  high  Mach  number  at  the  exit. 
This  suggests  that  cumulative  nonlinear  acoustic  distortions  are  unlikely 
to  occur  under  subsonic  jet  exit  conditions.  Albeit,  viscous  absorption 
is  considerably  weaker  at  low  frequencies,  and  so  it  is  still  plausible  for 
a  geometrically  large  subsonic  flow  to  produce  acoustic  waves  capable  of 
undergoing  cumulative  nonlinear  distortions.  Nevertheless,  Mach  waves 
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that  develop  in  supersonie  jet  flows  produce  pro-steepened  waveforms  at 
the  source  which  can  accelerate  the  wave  steepening  process  and  shorten 
the  shock  formation  distance.  A  carefully  documented  execution  of  the 
generalized  Burgers  equation  using  pre-steepened  waveforms  could  pro¬ 
vide  some  useful  answers  to  this  concern. 

■  Currently,  the  model  computes  the  shock  formation  distance  of  the  broad¬ 
band  noise  as  if  the  sound  pressure  along  the  peak  noise  path  obeys  by 
spherical  spreading  (p  oc  1/p).  However,  in  close  vicinity  to  the  emission- 
point  location,  the  sound  pressure  decays  cylindrically  (p  oc  1  /y/p)  due 
to  the  non-compact  nature  of  the  noise  sources.  And  so,  when  taking 
that  into  account  accordingly,  the  shock  formation  distance  will  increase. 
However,  a  pre-steepened  waveform  addressed  in  the  previous  point,  will 
result  in  a  shorter  shock  formation  distance. 

Temporal  waveforms  from  the  Mach  3  jet  were  examined  using  various 
statistical  (ensemble- aver  aged)  metrics  including  the  skewness,  kurtosis,  wave 
steepening  factor,  number  of  zero  crossings  and  the  Morfey-Howell  nonlinearity 
indicator.  Contours  of  these  metrics  in  the  ( x ,  r)-plane  reveal  that  the  process 
of  nonlinear  acoustic  wave  steepening  is  absent  within  the  measurement  range 
considered.  That,  in  combination  with  the  findings  from  the  aforementioned 
scaling  model,  provides  confidence  m  stating  that  the  detection  of  significant 
cumulative  nonlinear  waveform  distortions  in  common  laboratory-scale  (range- 
restricted)  environments  is  practically  impossible.  It  is  thus  precluded  that 
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observed  shock-structures  are  the  consequence  of  nonlinear  wave  steepening 
over  distance.  As  a  final  note,  it  is  advocated  that  when  the  Morfey-Howell 
indicator  is  used  as  an  investigative  tool  for  the  presence  of  cumulative  nonlin¬ 
ear  waveform  distortions,  it  must  be  applied  as  a  multi-point  indicator  along 
spatial  paths  along  which  the  sound  travels  in  an  ensemble-averaged  fashion. 
Only  when  the  statistical  nature  of  the  source  is  known  to  be  purely  Gaussian, 
the  indicator  can  be  applied  as  a  single-point  indicator.  However,  this  is  not 
applicable  to  jets,  since  the  statistical  nature  of  the  noise  sources  are  unknown. 

7.2  Crackle  emitted  by  high-speed  jets 

Several  methods  were  pursued  to  quantify  acoustic  shock-structures  in 
the  sound  field  of  the  experimental  Mach  3  jet.  The  absence  of  any  cumu¬ 
lative  nonlinear  steepening  effects  in  the  measurement  range  (§  7.1)  implied 
that  acoustic  shock-structures  are  generated  by  local  mechanisms  in,  or  in 
close  vicinity  to,  the  jet’s  hydrodynamic  region  and  that  they  can  be  classi¬ 
fied  as  crackle.  Some  believe  that  the  separation  between  the  hydrodynamic 
region  and  acoustic  region  of  a  jet  is  equivalent  to  the  separation  identified 
by  rotational  and  irrotational  flow.  The  research  on  crackle  was  motivated  by 
the  fact  that  (1)  it  is  perceived  as  the  most  annoying  component  of  jet  noise, 
(2)  no  unique  measures  of  crackle  exist,  and  (3)  significant  reductions  in  jet 
noise  will  be  achieved  when  crackle  can  be  controlled.  Considering  the  second 
point,  only  a  few  statistical  metrics  have  been  applied  to  pressure  waveforms 
in  the  past  with  an  attempt  to  assess  the  existence  of  crackle.  Shortcomings  of 
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using  these  ensemble- averaged  metrics  are  that  they  do  not  solely  weight  the 
waveform’s  shock  content,  blit  rather  take  the  total  waveform  into  considera¬ 
tion.  This  acts  as  a  factor  of  noise  in  the  quantification  process.  Henceforth, 
crackle  is  better  quantified  by  studying  the  shock  content. 

A  shock  detection  algorithm  (SDA)  was  developed  to  detect  shock 
fronts  in  the  acoustic  pressure  waveforms.  The  algorithm  was  based  on  two 
user  defined  thresholds  of  the  temporal  pressure  derivative  (pressure  rise  time) 
and  the  pressure  jump  across  the  shock  (shock  strength).  At  first,  the  average 
shapes  of  the  waveforms  associated  with  crackle  were  obtained  by  ensemble¬ 
averaging  all  instances  where  shock-structures  had  been  identified  using  the 
SDA.  Shocks  are  close  to  being  point-symmetric  around  the  shock  front  (de¬ 
fined  as  the  median  time  between  pressure  minima  and  maxima)  at  shallow 
angles  to  the  jet  axis  and  have  rounded  pressure  minima  and  maxima.  At 
angles  near  the  Mach  cone  edge,  the  pressure  maxima  become  more  sharp  and 
are  indicative  of  less  pronounced  shock  relaxation  effects.  It  was  concluded 
that  classifying  the  shock-structures  identified  outside  of  the  Mach  cone  as 
being  crackle  is  highly  questionable.  Furthermore,  within  the  Mach  cone,  the 
average  number  of  shocks  per  second  did  not  vary  along  radial  lines  of  constant 
polar  angle  0  with  a  focal  point  residing  in  the  region  of  most  intense  sound 
generation  (x  =  17.5  D^)  And  so  the  crackling  structures  travel  efficiently  in 
the  sense  that  they  do  not  disappear  due  to  relaxation  within  the  measurement 
range.  Trends  of  temporal  shock  instances  (arrival  times)  were  subsequently 
studied.  Like  the  average  number  of  shocks  per  second,  PDF  s  of  the  l5t-order 
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intermittence  are  independent  of  the  position  on  each  of  the  radial  lines  A- 
G.  The  revealed  modal  behavior  of  the  time  instances  provides  evidence  that 
shock-structures  appear  to  be  present  in  groups  of  multiple  shocks.  It  was 
postulated  that  physical  mechanisms  leading  to  these  modal  observations  are 
related  to  different  families  of  curved  wave  fronts  that  form  fringe  patterns. 
Suggestions  for  future  improvements  in  the  detection  and  quantification  of 
temporal  shocks  are  itemized  as  follows. 

■  Currently,  the  SDA  is  based  on  a  fixed  rise  time  (the  ratio  of  the  two  user 
defined  thresholds,  A ts  =  0.0540  ms).  However,  from  studies  by  Pierce 
[64]  that  are  reviewed  in  the  work  of  Mclnerny  &  Oilmen  (2005)  [46].  it 
became  apparent  that  for  stationary  weak  shocks  (ones  that  propagate 
without  a  change  in  shape),  the  rise  times  are  related  to  the  jumps  in 
pressure.  And  so,  the  SDA  shock  selection  criteria  can  be  adapted  to 
reflect  this  dependence. 

■  Wavelet- based  time- frequency  analyses  are  used  throughout  the  liter¬ 
ature  for  shock  detection  [66,  1].  In  this  research,  the  wavelet  power 
spectra  are  used  secondary  to  the  SDA  to  reveal  spectral  signatures  of 
crackle.  In  future  work,  wavelet-based  SDA  can  provide  extra  confidence 
in  the  correctness  of  shock  detection. 

■  Various  concepts  explaining  the  temporal  shock  intermit tencies  were  ad¬ 
dressed  in  a  preliminary  fashion  (§  6.3.2).  Synchronized  measurements 
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of  shock-structures  (acoustic  pressure)  and  pressure  or  velocity  measure¬ 
ments  in  the  hydrodynamic  region  of  the  jet  (preferably  spatially  exten¬ 
sive)  arc  believed  to  shed  more  light  on  these  ideas.  Furthermore,  fu¬ 
ture  work  with  analytical  models  encompassing  the  propagation  of  wave 
fronts  can  be  insightful  to  reveal  the  true  physical  mechanisms  of  how 
these  crackle-inducing  structures  manifest  themselves  in  the  sound  field. 
Wave  packet  modeling  (appendix  A)  appears  to  be  a  promising  tech¬ 
nique,  especially  when  a  sequence  of  wave  packets  is  superposed,  each 
having  a  unique  axial  wave  number  and  pressure  amplitude  envelope. 

A  spectral  measure  of  the  level  of  crackle  was  considered  by  using 
wavelet-based  time- frequency  analyses.  By  employing  the  wavelet  transform 
of  acoustic  waveforms,  acoustic  pressure  spectra  are  generated  as  function  of 
time.  Thereafter,  the  time- frequency  energy  distribution  is  used  to  detect 
the  unique  spectral  content  associated  with  the  shock  fronts  identified  by  the 
SDA,  this  is  denoted  as  the  shock  spectrum.  In  order  compute  a  practically 
significant  measure,  laboratory-scale  spectra  are  scaled  to  an  equivalent  full- 
scale  scenario  at  first.  Secondly,  spectra  are  A- weighted  to  account  for  the 
relative  loudness  perceived  by  the  human  ear  at  different  frequencies  within 
the  audible  range.  The  absolute  increase  of  spectral  sound  energy  between  the 
global  jet  noise  spectrum  and  the  shock  spectrum  is  independent  of  outward 
distance  along  radial  lines  and  complies  with  the  findings  of  constant  tempo¬ 
ral  characteristics  along  the  radial  paths  with  constant  polar  angle  cj).  The 
increase  in  sound  energy  is  quantified  as  the  percentage  energy  gain  from  the 
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global  pressure  spectrum  to  the  shock  spectrum.  This  energy-based  metric  is 
postulated  to  be  an  appropriate  metric  for  the  perception  of  crackle  since  it 
provides  the  correct  spatial  trends:  crackle  is  perceived  dominantly  at  shallow 
angles  and  gradually  becomes  less  pronounced  at  larger  angles.  Suggestions 
for  future  work  on  the  spectral  measures  of  crackle  are  itemized  below. 

■  In  order  to  conclude  that  energy-based  quantification  of  crackle  is  rele¬ 
vant  for  the  perception  of  crackle,  a  psycho-acoustic  study  involving  blind 
jury  trials  has  to  be  performed.  The  subjectiveness  of  crackle  plays  an 
important  role  in  a  quantitative  measures  of  ones  ability  to  reduce  jet 
noise. 

■  On  a  final  note,  crackle  and  cumulative  nonlinear  distortions  are  most 
likely  related.  Findings  suggest  that  crackle  is  not  the  consequence  of 
cumulative  nonlinear  waveform  distortions.  However,  cumulative  effects 
can  aid  in  the  effectiveness  of  crackle,  meaning  that  these  effects  pre¬ 
vent  early  relaxation  of  crackling  structures,  which  therefore  enhance 
the  amplitude  and  perception  of  crackle. 
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